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ABSTRACT 

We find that at redshifts z 10, HD line cooling allows strongly-shocked primordial 
gas to cool to the temperature of the cosmic microwave background (CMB). This 
temperature is the minimum value attainable via radiative cooling. Provided that the 
abundance of HD, normalized to the total number density, exceeds a critical level of 
~ 10 -8 , the CMB temperature floor is reached in a time which is short compared 
to the Hubble time. We estimate the characteristic masses of stars formed out of 
shocked primordial gas in the wake of the first supernovae, and resulting from the 
virialization of dark matter haloes during hierarchical structure formation to be ~ 
IOMq. In addition, we show that cooling by HD enables the primordial gas in relic H II 
regions to cool to temperatures considerably lower than those reached via H2 cooling 
alone. We confirm that HD cooling is unimportant in cases where the primordial gas 
does not go through an ionized phase, as in the formation process of the very first 
stars in z <; 20 minihaloes of mass ~ 1O 6 M . 
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1 INTRODUCTION 

What were the properties of the first generations of stars 
that formed at the end of the cosmic dark ages? Within the 
standard picture of hierarchical structure formation, shocks 
in the primordial gas arise naturally, both through the merg- 
ing of dark matter (DM) haloes in the formation of galaxies 
and through the supernovae (SNe) that marked the deaths 
of the first generation stars (e.g. Barkana & Loeb 2001; Cia- 
rdi & Ferrara 2005). During the merging of DM haloes, the 
primordial gas embedded in them can be shock-heated to 
high temperatures (<; 10 4 K), owing to the large velocities 
with which the haloes collide. Blast waves, triggered by the 
first SNe in the early Universe, are expected to have been 
commonplace (e.g. Madau, Ferrara & Rees 2001; Mori, Fer- 
rara & Madau 2002; Bromm, Yoshida & Hernquist 2003), 
as the first generation of stars were likely very massive 
(Bromm, Coppi & Larson 1999, 2002; Abel, Bryan & Nor- 
man 2002; Nakamura & Umemura 2001). The first stars, 
called Population III (Pop III) stars because they are as- 
sumed to have near-zero metallicity (e.g. Bromm & Lar- 
son 2004; Glover 2005), are thus predicted to be sufficiently 
shortlived (Si 3 Myr) to explode within a time short com- 
pared to the Hubble time (e.g. Bond, Arnett & Carr 1984). 
It has long been realized that the evolution of shocked 
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primordial gas is markedly different from the nearly adia- 
batic one in the early minihaloes with total masses ~ 1O 6 M0 
and collapse redshifts of z <; 20 inside of which the very 
first stars are predicted to have formed (e.g. Mac Low & 
Shull 1986; Shapiro & Kang 1987; Yamada & Nishi 1998; Oh 
& Haiman 2002; Clarke & Bromm 2003; Omukai & Yoshii 
2003). In the strong-shock limit, the shocked, ionized gas 
cools roughly isobarically, and can cool faster than it recom- 
bines. Under these conditions, molecular hydrogen is readily 
formed, and the resultant high abundance of H2 allows the 
gas to cool to Si 200 K and to achieve high densities (e.g. 
Kang & Shapiro 1992). The cooling of primordial gas be- 
hind shock fronts thus allows gas fragmentation to occur 
at lower temperatures, possibly leading to the formation of 
stars with masses that are smaller than in the unshocked 
case (e.g. Mackey, Bromm & Hernquist 2003; Salvaterra, 
Ferrara & Schneider 2004; Machida et al. 2005). 

The presence of deuterium in the primodial gas has re- 
cently received renewed attention as a means to allow the 
gas to cool to temperatures yet lower than those possible 
solely via H2 line cooling (e.g. Flower et al. 2000; Galli & 
Palla 2002; Nakamura & Umemura 2002; see also Suchkov, 
Shchekinov & Edelman 1983). The role of HD in star forma- 
tion, through the fragmentation of SN-shocked primordial 
gas, has been studied as well, suggesting that line cooling 
by HD may facilitate the formation of low mass metal-free 
stars, possibly even of primordial brown dwarfs (Uehara & 
Inutsuka 2000; Machida et al. 2005). 
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very small or zero 
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£ 10 8 M o haloes 


II 


~ 10- 2 


~ 1 


Early galaxies 


I 


£0.1 


~ 1 


Present-day galaxies 



Table 1. Stellar populations in the high-redshift Universe. For each population, we show the typical mctallicity of the parent gas, the 
estimated fragmentation mass scale, and the typical formation site (see text for a detailed discussion). 



Employing the extensive deuterium chemical network 
compiled by Nakamura & Umemura (2002), we investigate 
the importance of HD in the cooling of primordial gas, fo- 
cusing on four distinct cases: shocks from SNe, shocks as- 
sociated with the virialization of DM haloes in large-scale 
structure formation, collapse within minihaloes, and collapse 
within relic H II regions. In addition, we provide estimates 
for the characteristic mass of stars formed by the fragmenta- 
tion of strongly-shocked primordial gas. The outline of our 
paper is as follows. In Section 2 we discuss the basic phys- 
ical ingredients included in our treatment of the cooling of 
shocked primordial gas. Section 3 describes the evolution 
of the primordial gas in structure formation shocks and in 
SN-induced shocks, as well as in minihaloes and relic H II 
regions. In Section 4 we derive the characteristic mass of 
stars formed through the fragmentation of shocked primor- 
dial gas. Finally, we present our conclusions in Section 5. 



2 BASIC PHYSICAL INGREDIENTS 

2.1 Shocks in the high-z Universe 

We consider two different situations in the early universe 
where the primordial gas is strongly shocked: the assembly 
of DM-dominated haloes in large-scale structure formation 
and SN explosions of the first stars. In each of these cases 
the primordial gas is shock-heated to high temperatures and 
ionized, after which it cools roughly isobarically (Shapiro & 
Kang 1987; Yamada & Nishi 1998). This condition allows 
the gas to become very dense as it cools behind the shock 
front, as well as allowing it to cool more quickly than it can 
recombine. Thus, the gas remains out of equilibrium even 
after cooling to below ~ 10 4 K, and the formation of molec- 
ular hydrogen is enhanced, catalyzed by free electrons and 
hydrogen ions (see Kang & Shapiro 1992). The result is an 
enhanced fraction of H2, the line cooling of which can then 
efficiently lower the gas temperature all the way down to 
<, 200 K. Without sufficient H2, gas of primordial composi- 
tion can obtain, via atomic hydrogen line cooling, tempera- 
tures no lower than ~ 10 a ' 9 K. 

The high densities and low temperatures attainable 
by the shocked gas in this way, through isobaric cooling 
and the effective out-of-equilibrium formation of molecu- 
lar hydrogen, make possible a much lower Jeans mass than 
would otherwise be the case for a primordial gas cloud, since 
Mj oc T 3 / 2 n _1,/2 , where T is the temperature of the gas and 
n is the number density of its constituent particles. This, in 
turn, translates into lower characteristic masses for the stars 
that form from such shocked primordial gas. 

This is one mechanism that would lead to the formation 
of so-called Population II. 5 (Pop II. 5) stars, which have been 



postulated to have formed from the primordial gas shocked 
by the SN explosions of Pop III stars (Mackey et al. 2003; 
Salvaterra et al. 2004). More generally, Pop II. 5 stars could 
have formed from primordial gas that has been strongly ion- 
ized, thus allowing for the free electron-catalyzed formation 
of molecules and so for enhanced cooling of the primordial 
gas, regardless of the source of the ionization. Another pos- 
sible site for the formation of Pop II. 5 stars thus are DM 
haloes of masses i> 10 8 M©, containing gas shocked through 
the hierarchical merging of DM haloes in the process of early 
structure formation, at redshifts z J> 10. Finally, Pop II. 5 
stars may also have formed in the relic H II regions of the 
first stars, as the formation of molecules would have been 
catalyzed by the residual fraction of free electrons in the 
cooling gas following the death of the central massive star 
(e.g. O'Shea et al. 2005). Pop II. 5 stars would have been of 
almost purely primordial composition and with masses at 
least an order of magnitude lower than the typical masses of 
Pop III stars. A careful treatment of the primordial chem- 
istry, including the cooling effects of HD, appears to be nec- 
essary in determining the characteristic masses of Pop II. 5 
stars. 

To summarize, the postulated Pop II. 5 stars, defined 
as stars formed from strongly ionized primordial gas, have 
masses intermediate between Pop III and Pop II ones, but 
may still be almost metal-free. The distinct modes of star 
formation are discussed in Table 1. 

2.2 Chemistry of the primordial gas 

In evaluating the cooling properties of shocked primordial 
gas, we adapt the one-zone model used in Mackey et al. 
(2003). This model solves for the thermal and chemical evo- 
lution of isobarically cooling gas of primordial composition, 
though insofar only as including normal hydrogen and he- 
lium (see Mackey et al. 2003; Bromm et al. 2002). In ad- 
dition to these chemical species, we have included for our 
calculations the five deuterium species D, D + , D~, HD, and 
HD + , together with a network of 18 reactions in which these 
species participate. The rate coefficients for these reactions 
are given in table 1 of Nakamura & Umemura (2002). We 
replace the rate coefficients for the following reactions with 
the updated rate coefficients from Galli & Palla (2002): 

D + H+ -f D+ + H , 

D+ + H -f D + H+ , 

D+ + H 2 ->■ H+ + HD , 

HD + H + -> H 2 + D+ . 

We would like to point out that the rate coefficient for 
reaction (3), as quoted in Galli & Palla (2002), contains a 
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Figure 1. The equilibrium abundances of the major species in- 
cluded in calculating the evolution of primordial gas. The solid 
lines show the equilibrium abundances to which our chemical net- 
work converged, while the diamond points are the results from 
analytical calculations of the equilibrium abundances of the two 
most important coolants of the gas at low temperatures, H2 and 
HD. The electron abundance exceeds unity at high temperatures 
because of the contribution of electrons due to the ionization of 
helium at those high temperatures. The dashed line in the panel 
for the electron abundance denotes the maximum abundance pos- 
sible for the electrons, which occurs only when helium is com- 
pletely ionized. 

typo, and we refer the reader to the original reference for 
this rate (Savin 2002) for the correct value. 

We expect that the most important coolant in the pri- 
mordial gas at temperatures ^ 100 K is HD, and hence we 
neglect the presence of other trace species which could have 
formed in the primordial gas, such as LiH owing to the very 
low primordial abundance of lithium (n^i/n ~ 10~ 10 ). We 
also neglect the cooling of the primordial gas due to ionic 
molecules such as H2 + , H3 + , and HeH + , which could be 
significant in the temperature range 5,000 to 8,000 K. Ex- 
cluding these coolants only strengthens our result, since any 
additional cooling at < 8, 000 K would enhance the out-of- 
equilibrium production of H2 and HD. 

Our chemical network was tested by verifying that the 
equilibrium abundances of the main coolant species, as cal- 
culated numerically, matched those derived analytically. The 
results of these calculations are displayed in Fig. 1, the solid 
lines showing the equilibrium abundances converged upon 
by the code and the diamonds representing the analytical 
ones. We limit our analytical derivations to the two most 
important coolants at the low temperatures at which the 
gas fragments 200 K), H2 and HD. As Fig. 1 shows, the 
numerical results nicely reproduce the analytical values (for 
details of the analytical estimates see Bromm et al. 2002) . 

The primary reaction sequence which creates HD in the 
shocked primordial gas is (e.g. Galli & Palla 2002; Machida 
ct al. 2005) 

e" + H^H~ +hv , (1) 
H-+H^H2+e", (2) 




10 2 10 3 10 4 10 5 

Temperature [K] 



Figure 2. The abundances of HD, H2, D, and D + as a function 
of temperature, corresponding to the 100 km s — 1 shock at z = 
20 (see lower left panel of Fig. 3). The dominant reaction forming 
HD is D+ + H 2 — > HD + H+. Despite the steady decline in the 
abundance of D+ with decreasing temperature, HD is still created 
predominantly through the above reaction. 

D+ + H 2 -> HD + H+ . (3) 

The dependence of the HD abundance on those of H2 and 
D + is illustrated in Fig. 2, which shows the abundances of 
H2, HD, D, and D + as a function of temperature in the 
case of a SN-induced shock, corresponding to the lower left 
panel of Fig. 3. Evident from Fig. 2 is a steady increase in 
HD abundance with lower temperature, and so with higher 
density, at temperatures T gas 500 K, which is due to the 
out-of-equilibrium abundance of H2 and the residual abun- 
dance of D + which has not recombined even at the lowest 
temperatures T gas ~ Tcmb • 

2.3 Photodestruction of HD 

It is well established that photons in the Lyman- Werner 
(LW) bands, with energies between 11.2 and 13.6 eV, can 
easily dissociate H2 molecules (e.g. Haiman, Rees & Loeb 
1997; Bromm & Larson 2004 and references therein). Here 
we evaluate the effect that such radiation has in dissociating 
HD molecules. 

Regardless of the source of the radiation, the timescale 
for the photodissociation of HD by LW photons can be ap- 
proximated by fdiss — 1000/x yr, where \ is the LW radia- 
tion field in units of the Draine standard field (see figure 1 
in Petit, Roueff & Le Bourlot 2002). We obtain the value 
of the Draine standard LW field by integrating the photon 
flux defined by equation (11) of Draine (1978) from 11.2 eV 
to 13.6 eV, which yields ~ 10 -20 erg s _1 cm -2 Hz -1 sr _1 . 
Thus, for a general radiation field the HD photodissociation 
time can be expressed as idiss ~ 10 8 yr (Jlw/IO -4 ) -1 , where 
Jlw measures the LW flux in units of 10 -21 erg s _1 cm~ 2 
Hz" 1 sr- 1 . 

To estimate the importance of radiation emitted by 
shocked primordial gas in destroying HD molecules, we com- 
pare the HD photodissociation timescale, tdiss, to the for- 
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mation timescale, if orm (see Machida et al. 2005). The LW 
photon flux emergent from a 100 km s~ shock passing 
through a medium of pre-shock hydrogen number density 
of 10 2 cm~ 3 is ~ 10 2 photons cm~ 2 s _1 , which corresponds 
to J LW ~ 10~ 4 (Shull & Silk 1979). This results in t diss 
~ 10 8 yr (Le Petit, Roueff & Le Bourlot 2002). From our 
calculations, we estimate t{ OIln ~ 10 6 yr for the same shock 
velocity and pre-shock density. Thus, for cases in which the 
post-shock gas is the only source of LW photons, the forma- 
tion timescale is shorter than the dissociation timescale by 
two orders of magnitude, and the effect of photodissociating 
radiation produced by the shock is likely to be unimportant 
over the range of shock velocities and pre-shock densities 
that we consider in this paper. 

We note, however, that in relic H II regions the pho- 
todestruction of HD may be particularly important due to 
emission of LW radiation through the two-photon decay of 
recombining helium atoms, at densities below ~ 10 4 cm~ 3 
(Mathis 1957; Pottasch 1961). While helium recombination 
radiation will also be generated by shocks strong enough to 
ionize helium, the column densities of ionized helium in the 
first relic H II regions are expected to be much higher than 
those in shock fronts, since the first H II regions generated 
by massive Pop III stars are expected to have been of order 
a few kpc in diameter with the majority of the helium be- 
ing ionized by the hard UV emission from the central star 
(Kitayama et al. 2004; Whalen, Abel & Norman 2004; Al- 
varez, Bromm & Shapiro 2005). With these larger column 
densities, the fluxes of recombination radiation from helium 
in relic H II regions may be strong enough to efficiently pho- 
todissociate HD molelcules. However, if molecule formation 
takes place efficiently after the gas has cooled below the 
temperatures at which helium recombines, then the effect of 
the destruction of molecules by helium recombination radia- 
tion could be mitigated. A careful treatment of the effect of 
recombination radiation is clearly required in order to accu- 
rately calculate the chemical and thermal evolution of relic 
H II regions. We defer such a treatment to future work. 

If the column density of HD is high enough that self- 
shielding is important, the value of t diBB calculated here will 
represent a conservative lower limit on the photodissocia- 
tion time of shielded HD molecules. The effects of H2 self- 
shielding have been widely studied and have been found to 
have important consequences for the cooling of primordial 
gas subject to a LW background radiation field (Draine & 
Bertoldi 1996; Haiman et al. 1997; Yoshida et al. 2003). As 
well, it has been shown that the dissociation rate of HD 
shielded by H2 is expected to be lower than the dissocia- 
tion rate of unshieldded HD by up to a factor of about 1/3, 
due to the close overlapping of several HD lines by H2 lines 
(Barsuhn 1977). Here, for simplicity, we neglect the effects 
of both HD self-shielding and HD shielding by H2. We note, 
however, that if there exists a general LW background in ad- 
dition to the radiation emitted locally by the shocked gas, 
such as from first generation massive stars during structure 
formation or from continuum radiation from recombining 
helium atoms in relic H II regions, then self-shielding ef- 
fects may become much more important in preventing the 
photodissociation of HD. 



2.4 Cooling by HD 

The cooling properties of HD are distinct from those of H2 
for the following reasons (e.g. Flower et al. 2000; Uehara & 
Inutsuka 2000). Firstly, HD has a permanent dipole moment, 
allowing dipole rotational transitions, which spontaneously 
occur much more often than the quadrupole rotational tran- 
sitions in H2. Also, the dipole moment of HD allows rota- 
tional transitions of AJ = ±1, which are of lower energy 
than the A J = ±2 quadrupole transitions of H2. Thus, col- 
lisional excitation of HD from the ground to the first excited 
rotational level (J=l), and the subsequent radiative decay 
back to the ground state by a dipole transition, can allow 
HD to efficiently cool gas to lower temperatures than can be 
reached with H2 line cooling alone. Finally, the rotational 
energy levels of HD are lower than those of H2 by a factor 
of the reduced mass, which is higher by a factor of ~ 4/3 for 
HD than for H2, allowing for even lower energy collisions to 
excite the first rotational level. 

We have implemented the HD cooling function provided 
by Flower et al. (2000; see also Lipovka et al. 2005). To 
take into account that the gas cannot radiatively cool to 
below the temperature of the cosmic microwave background 
(CMB), 



Tcmb = 2.7 K(l + z) 



(4) 



we employ an effective cooling rate, according to 

A = A(Tga B ) — A(Tcmb) , (5) 

where T gas is the temperature of the gas. 

2.5 The CMB cooling limit 

For a gas which cools purely radiatively, and not by adia- 
batic expansion, the CMB sets a lower limit on the tem- 
perature to which it can cool (e.g. Larson 1998). To show 
this, we will consider the process by which the primordial 
gas cools to temperatures approaching the CMB value. The 
frequency of emitted radiation for the rotational transition 
J = 1 — » of HD is given by 



hvio 
kB 



130 K 



(6) 



where vio is the frequency for the J = 1 — * rotational 
transition, fee is the Boltzmann constant, and h is the Planck 
constant. We investigate how the primordial gas approaches 
Tcmb by only considering this transition and its reverse. 
This is reasonable since all other transitions from the ground 
state of HD have excitation energies significantly higher than 
that for this transition and since, for the redshifts we are 
considering, Tcmb < 130 K. 

Denote the Einstein coefficients for spontaneous and 
stimulated emission from J = 1 — > by A w and Bio, re- 
spectively. The Einstein coefficient for absorption of a pho- 
ton effecting the transition J = — > 1 is B01 • Now, consider 
a finite parcel of optically thin primordial gas at a tempera- 
ture T gas . We assume for simplicity that densities are suffi- 
ciently high to establish LTE level populations according to 
the Boltzmann distribution 



m 
no 



9o 



(7) 
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where gi and go are the statistical weights of the first excited 
and ground rotational levels with g\ = 3<;o- Furthermore, it 
is assumed that only the ground and first excited rotational 
levels are occupied. That is, we take it here that already 
T gas < f30 K. We also make the assumption that T gas <; 
Tomb- Thus, if we denote the blackbody specific intensity of 
the CMB at the frequency vio as I vw , then it follows that 



Ihvlolc 2 



e fc B J CMB 



(8) 



g B ga 



If we use the standard relations Bioffi = Boigo and 

2hu'j _ Aio 
c 2 ~ Bio ' 

along with equation (7), then equation (8) can be written as 



(9) 



7 Aio / -going _ x > 
" 10 Bio V-Biom 



(10) 



For HD, the Einstein coefficient for spontaneous emission 
is A w - 5x 10~ 8 s _1 (e.g. Nakamura & Umemura 2002). 
This, in turn, implies that 



no 



BoiI vla < niAio + niB\oI vw 



(11) 



Thus, more energy is emitted into the radiation field 
than is absorbed from the radiation field. If we define the 
energy density of the gas according to 



^gas — ~n/cBT ga s , 



(12) 



where n is the total number density of the gas particles, 
including all species, then equation (11) implies that, with 
no change in the density of the gas, 

3 dT 

hv w [noBoiI vlo - niA w - mB 10 I vw ] = -nk B — ^ . (13) 

2 at 

Despite our assumption of LTE level populations, no 
collisional excitation and de-excitation terms appear in 
equation (13), because it expresses the exchange of energy 
between two physically distinct heat reservoirs, the CMB 
and our parcel of gas. In LTE, collisional excitations are ex- 
actly balanced by de-excitations, and there would be no net 
heat transfer between the two reservoirs due to these terms. 
Next, we take it that the ratio of the number density of HD 
molecules «hd to the total number density of particles n in 
the gas is given by the constant factor 



nun 
n 



no + m 



no 
n 



We note that, since Tomb Ss Tgas < 130 K, 

e *B T CMB <J e fc B T gas < 1 . 



This implies that 



2hv( 



lfl_ 



lfl_ 



Aio - - —r 

B J CMB = e 'B T CMB 

Bio 



(14) 



(15) 



(16) 



For simplicity, we further assume that stimulated emission 
can be neglected. The thermal evolution of the gas is then 
approximately described by 



dT ei 



dt 



Ihvio Aio-Xhd 

fcB 



ft "10 ""10 
fe B r CMB — e fc B T gas 



(17) 



It is clear from this result that if Tcmb ^ T gas < 130 K, 
with the gas cooling only by radiative decay of the excited 
rotational level J—l to J=0, the temperature of the gas 
will asymptotically approach Tcmb - Thus, equation (17) de- 
scribes the fact that the CMB temperature is indeed a lower 
limit on the temperature to which a gas can cool via line 
emission only. 

We estimate the timescale for reaching the CMB tem- 
perature floor to be 



£cmb ~ 1 ( feTcMB \ 2 ex / hvio \ 

2Ai X H d V hvio ) VfesTcMB/ 
— (^4io^hd) 1 • 

Defining a critical HD abundance by demanding tcMB 
with tu being the Hubble time, we find 



10" 



(18) 



(19) 



This estimate, derived for the limiting case of LTE level pop- 
ulations, is almost independent of redshift for z ij 10. There- 
fore, if HD has a fractional abundance lower than XHD.crit 
~ 10~ 8 , the gas will not cool to T gas ~ Tcmb within a Hub- 
ble time at redshifts z 10. Thus, a mechanism for produc- 
ing large fractional abundances of HD is necessary in order 
for the primordial gas to cool to Tcmb before fragmentation 



3 EVOLUTION OF PRIMORDIAL GAS 

3.1 First generation SN blast waves 

The passage of a SN shock through the primordial gas re- 
sults in the almost complete ionization of the hydrogen for 
shock velocities J> 100 km s _1 , early in the radiative phase 
of the shocks (see Kang & Shapiro 1992). Fig. 3 shows the 
thermal evolution of primordial gas compressed, heated, and 
ionized by shocks with velocities in the range expected for 
the first radiative supernova shocks, 100 and 200 km s~ , 
and at redshifts of z = 10 and z — 20. The initial post- 
shock temperature of the gas is defined by 



Tr>; 



3&B 



(20) 



where M sh is the velocity of the shock, and mn is the mass 
of a hydrogen atom. The post-shock gas initially cools by 
brcmsstrahlung, and then by He line emission, once the tem- 
perature has dropped enough for helium atoms to recom- 
bine. Once hydrogen begins to recombine the temperature 
of the gas can drop to ~10 4 K by atomic hydrogen line emis- 
sion. 

For the cases shown in Fig. 3, large fractions of H2 
and HD form out-of-equilibrium as the gas cools below 
~10 4 K, allowing the gas to efficiently cool to temperatures 
Tg as & 200 K. In each case shown in Fig. 3, the gas con- 
taining deuterium is then able to cool to a temperature T gas 
— Tcmb in S3 10 7 yr, a fraction of the age of the Universe 
at redshifts z ^ 10. However, without deuterium, cooling 
only takes place by H2 line emission at low temperatures (S> 
200 K), and so it takes at least ~10 2 times longer for the 
temperature of the deuterium-free gas to approach the CMB 
temperature. Thus, while HD can rapidly cool the primor- 
dial gas to the CMB temperature, H2 does not efficiently 
cool below <; 200K. 



© 2005 RAS, MNRAS 000, 000-000 



6 J.L. Johnson and V. Bromm 



z=10 
u a =100 km s"' 



-8 1Q -6 


10 4 10 2 10° 
Hubble times 




z=20 
u„=10O km s"' . 





z=10 
=200 km s" 1 



10 ° 10 ° 10 ' 10 ' 10" 
Hubble times 



z=20 

— s u„=200 km s" 1 




30 20 

Redshift [z] 



Figure 3. The thermal evolution of SN-shockcd primordial gas. 
The dashed lines show the temperature as a function of time for 
primordial gas without deuterium, while the solid lines show the 
gas temperature when deuterium is included. The combination 
of redshift z and shock velocity u s ^ is varied in each panel. In 
each case the temperature of the gas reaches the temperature 
of the CMB, denoted by the dotted line, in a small fraction of 
the Hubble time. Each case was calculated assuming a preshock 
density of no = 10 2 cm -3 . 



Although the thermal evolution in each of the cases 
plotted in Fig. 3 is calculated assuming a pre-shock density 
of no = 10 2 cm -3 , we have found that the gas containing 
deuterium is efficiently cooled to the CMB limit over a range 
of pre-shock densities. After the passage of a shock with a ve- 
locity of 200 km s _1 at z = 20, the gas cools to T gas c± Tcmb 
in ~ 10 _1 te for n pa = 10~ 2 cm~ 3 and in ~ 10~ 2 £h for no = 
1 cm~ 3 . Thus, for SN shocks propagating through primor- 
dial gas in a variety of situations, sufficient HD can form to 
allow cooling to the temperature of the CMB in a fraction 
of the age of the Universe at redshifts z J> 10. As shown in 
Fig. 2, the fractional abundance of HD at temperatures <, 
200 K is Xhd ii 10~ 6 , well above the minimum of XHD.crit 
~ 10 required for the effective cooling of the primordial 
gas by HD. 



3.2 High-2 structure formation shocks 

We here distinguish between the evolution of shocks associ- 
ated with SNe and with the assembly of DM halos during 
hierarchical structure formation. For the case of the cooling 
of primordial gas shock-heated during DM halo assembly, 
the virial velocities can be up to an order of magnitude lower 
than for SN shock velocities, while the resulting shock can 
still partially ionize the primordial hydrogen (see Shapiro & 
Kang 1987). 

Fig. 4 shows the virial velocities associated with the DM 
halo assembly for various halo masses. Using these veloci- 
ties and expressing the masses of the resultant DM haloes 
in terms of the a fluctuations they represent at a given red- 
shift (e.g. Barkana & Loeb 2001), Fig. 5 shows the ther- 
mal evolution of primordial gas shocked in the formation of 
haloes with masses between ~ 10 s and ~ 10 10 Mq at red- 



Figure 4. The virial velocities associated with the assembly of 
DM haloes during hierarchical structure formation. Each of the 
curves corresponds to a different halo mass. Because the densities 
of the virializcd haloes arc higher at higher redshifts, the virial 
velocities increase with redshift. 



shifts z = 10 and z — 15. We again assume that the gas 
is strongly-shocked, and so that the post-shock evolution 
is isobaric. The initial post-shock temperature is given by 
equation (20), and the pre-shock density of the gas is taken 
to be the density of baryons in DM haloes at the point of 
virialization, given by (e.g. Clarke & Bromm 2003) 



no ~ 0.3cm 



T + z 
. 21 



(21) 



If there are clumps of higher densities within the haloes, 
however, this estimate of the density will likely be too low. 
Then, because the gas can cool faster at higher densities, 
due to the more rapid formation of H2 and HD molecules, 
the cooling times shown in Fig. 5 may be overestimates for 
high density clumps within haloes. 

The initial post-shock fractional ionization of the var- 
ious chemical species composing the gas were taken from 
Shapiro & Kang (1987) and Kang & Shapiro (1992), for the 
cases of shock velocities of 20, 30, 50, and 100 km s _1 . As 
deuterium was not treated in this earlier work, we take it 
that the fractional abundances of D, D + , and HD are re- 
duced from those of H, H + , and H2, respectively, by a factor 
of the cosmic abundance of deuterium, taken to be ~ 4 x 
10~ 5 . 

As was demonstrated by Shapiro & Kang (1987), even 
at shock velocities 100 km s _1 , which characterize the 
collapse of DM haloes at redshifts 2 ^ 10, enough out- 
of-equilibrium formation of molecular hydrogen takes place 
that the gas can cool to 100 K within the age of the Uni- 
verse, at high enough redshift (see also Murray & Lin 1990). 
That H2 efficiently cools the gas behind these shocks is also 
shown in Fig. 5, where the temperature of the primordial 
gas does indeed fall to ~ 100 K within a Hubble time. The 
abundance of HD at low temperatures (jj 100 K) for the 
case of shocks due to DM halo collapse is found to be of the 
same order as that found for the SN case, that is Xhd iJ 
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Figure 5. The thermal evolution of primordial gas shocked 
through dark matter halo collapse. The dashed lines show the 
temperature as a function of time for primordial gas without deu- 
terium, while the solid lines show the temperature of the gas when 
deuterium is included in the calculation. The combination of red- 
shift z and a fluctuation of the DM halo is different in each panel. 
In each case the temperature of the gas reaches the temperature 
of the CMB, denoted by the dotted line, though generally not as 
quickly as is the case with the higher velocity supernova shocks. 
The initial preshock density at a given redshift was taken to be 
the virial density of baryons at that redshift. 

1CP 6 , again well above the minimum necessary for efficient 
HD cooling. 



Figure 6. The thermal evolution of an unshocked, collapsing pri- 
mordial gas cloud, such as would have resided in the first mini- 
haloes. Initially, the gas collapses adiabatically, starting with the 
density of baryons in a virializcd dark matter halo at a redshift 
of z = 20. Upon reaching a density of ~ 10 cm -3 , the gas cools 
by molecular rotational line-emission. The dashed line shows the 
temperature evolution of the primordial gas when deuterium is 
not included, leaving the gas to cool only by H2 line-cooling, 
which is effective at cooling the gas only down to ~ 200 K. The 
solid line shows the temperature evolution of the gas when deu- 
terium is included in the calculation. Although, in principle, HD 
can efficiently cool the gas to the temperature of the CMB, de- 
noted by the dotted line, unshocked gas of primordial composition 
does not contain sufficient HD (Xjjd ~ 10 -9 ) for it to be an ef- 
ficient coolant. 



3.3 Unshocked primordial gas within minihaloes 

Minihaloes, with masses of the order of 10 6 Mq, merge at 
velocities too low to strongly shock or ionize the primordial 
gas they contain, as can be seen from Fig. 4. Haloes of such 
low mass collide at velocities of only a few kilometers per 
second, which are not high enough to ionize the hydrogen 
inside the haloes, and so are not high enough to allow signif- 
icant free electron-catalyzed formation of molecules to take 
place. Thus, the gas in minihaloes is expected to remain 
largely un-ionized and with only the primordial fractions of 
H2 and HD until the formation of the first stars within them. 
Without the large fractions of HD formed in shocked, ionized 
primordial gas, the gas from which the first stars formed in 
minihaloes at redshifts of z Jj 20 could not cool below ~ 200 
K, and thus the first stars were likely very massive (Bromm 
et al. 1999, 2002). 

Fig. 6 shows the thermal evolution of unshocked gas of 
primordial composition with an initial ionization fraction of 
x c = 10~ 4 which collapses under its own gravity, its density 
evolving according to dn/dt = n/tg, where n is the number 
density of the gas and tg is the free-fall time, given by 

where p — m-an. Here, the initial abundances for this cal- 
culation were taken to be the primordial abundances of the 
standard cosmological model given by Galli & Palla (1998), 



at redshifts of z ~ 20. We take the initial temperature of the 
gas to be 200 K (e.g. Bromm et al. 2002) , and the initial den- 
sity of the gas is taken to be the density of baryons within 
DM haloes at the point of virialization, given by equation 
(21). 

The gas at first is heated by adiabatic compression as 
it collapses. Once reaching a density of ~ 10 cm~ 3 and a 
temperature of ~ 10 3 K, the gas begins to cool via H2 line 
emission, and the temperature drops again to ~ 200 K. As 
Fig. 6 shows, the fraction of HD stays too low for HD to 
efficiently cool the gas below these temperatures. We find 
that for almost a Hubble time after the collapse has begun, 
the abundance of HD has a value of Xhd 10~ s , which is 
below the minimum abundance necessary for effective HD 
cooling. Because the gas in this case is never ionized, the 
reactions (l)-(3), through which HD forms in the case of 
strongly-shocked primordial gas, are suppressed and HD is 
not effectively produced. Thus, cooling proceeds largely only 
through H2 line emission and the gas temperature never 
drops below ~ 200 K. 

3.4 HD cooling in relic H II regions 

The first stars which formed in minihaloes were likely 
strong sources of ionizing radiation due to their high mass 
(e.g. Bromm & Larson 2004). The possibility of a second- 
generation of stars forming in the so-called relic H II regions, 
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Figure 7. The thermal evolution of a primordial relic H II re- 
gion, in which an initially ionized primordial gas cloud collapses 
under its own gravity. The initial temperature of the gas was 
set to T gas = 3 X 10 4 K, typical for an H II region surrounding 
a massive Pop III star at z ~ 20. The dashed line denotes the 
temperature of the gas as a function of its increasing density, ne- 
glecting HD. The solid line shows the temperature evolution of 
the gas with HD cooling included. The free electrons, left to re- 
combinc following the extinguishing of the ionizing radiation from 
the hot star, catalyze the formation of H2, and consequently of 
HD. The abundance of HD becomes high enough to cool the gas 
to close to the CMB temperature at z = 20. 



carved out by the copious amounts of ionizing radiation 
emitted by first generation massive stars, has recently re- 
ceived some attention. In particular, it has been argued that 
the ionized conditions in relic H II regions allow for out-of- 
equilibrium production of large fractions of H2 (O'Shea et al. 
2005). Given that out-of-equilibrium fractions of HD can be 
created from H2 through reactions (l)-(3), large fractional 
abundances of HD may be created in relic H II regions, as 
well. Thus, HD may be an important coolant of the primor- 
dial gas in such regions. Fig. 7 shows the thermal evolution 
of primordial gas in a relic H II region. The initial temper- 
ature of the gas is taken to be 3 x 10 4 K, a temperature 
typical of an H II region surrounding a massive Pop III star 
(e.g. Storey & Hummer 1995). We thus take the hydrogen 
in the gas to be initially fully ionized and the fractional 
abundances of H2 and HD to be initially zero. The initial 
density, as in the cases of gas shocked in DM halo assem- 
bly and gas within minihaloes, is taken to be the density of 
baryons within DM haloes at the point of virialization, given 
by equation (21), for a redshift of z ~ 20. This density is 
consistent with the density of <, 1 cm~ 3 expected to result 
from the dynamical expansion of the first H II regions due 
to photoheating from a central massive star (see Kitayama 
et al. 2004; Alvarez et al. 2005). As in the case of unshocked 
gas within minihaloes, the evolution of the density n of the 
gas is described by dn/dt — n/ts- 

The effect of the out-of-equilibrium free electron- 
catalyzed production of H2 is evident, as the temperature of 
the primordial gas lacking deuterium falls to slightly lower 
values than in the case of unshocked, un-ionized primordial 



Figure 8. The evolution of the abundance of HD, Xhd> in pri- 
mordial gas which cools in four distinct situations. The solid line 
corresponds to gas with an initial density of 100 cm~ 3 which is 
compressed and heated by a supernova shock with velocity « s h 
= 100 km s^ 1 at z = 20. The dotted line corresponds to gas at 
an initial density of 0.1 cm~ 3 shocked in the formation of a 3cr 
fluctuation dark matter halo at z = 15. The dashed line corre- 
sponds to unshocked, un-ionized primordial gas with an initial 
density of 0.3 cm -3 collapsing inside a minihalo at z = 20. Fi- 
nally, the dash-dotted line shows the fraction of HD in primordial 
gas collapsing from an initial density of 0.3 cm -3 inside a relic H 
II region at z = 20. The temperature evolution in these situations 
is illustrated in figures 3, 5, 6, and 7, respectively. The horizontal 
line at the top denotes the cosmic abundance of deuterium. The 
critical abundance of HD, X[iD,crit> denoted by the bold dashed 
line, is that above which primordial gas can cool to the CMB 
temperature within a Hubble time, as given by equation (19). Al- 
though the HD abundance in the relic H II region does exceed 
Xhd criti adiabatic heating prevents the temperature of the gas 
in the relic H II region from cooling to the CMB temperature. The 
supernova-shocked gas and the gas shocked in the merging of DM 
haloes evolve isobarically and so adiabatic heating has much less 
effect on the temperature evolution in these situations. 

gas shown in Fig. 6. The larger fraction of H2 in the relic H 
II region case results in more effective H2 cooling. In turn, 
a large out-of-equilibrium abundance of HD results, as well, 
with the fractional abundance of HD exceeding Xhd ~ 10~ 5 
at the lowest temperatures (5s 100 K). With this enhanced 
fraction of HD the gas is able to cool almost to the temper- 
ature of the CMB at z ~ 20, despite the adiabatic compres- 
sion that heats the gas at high densities. 

Fig. 8 summarizes the evolution of the fractional abun- 
dance of HD, Xhd, for the cases previously described. HD 
cooling is significant whenever Xhd > -^HD.crit is reached in 
a fraction of the Hubble time. 

4 FRAGMENTATION OF PRIMORDIAL GAS 

4.1 Characteristic stellar mass for SN-shocked gas 

The characteristic mass of Pop II. 5 stars formed from SN- 
shocked primordial gas can now be estimated, assuming that 
the shock is strong, so that the shocked gas evolves isobari- 
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cally, and that the gas cools to the temperature of the CMB. 
Following Mackey et al. (2003), we assume that the imme- 
diate progenitor of a star is a self-gravitating gas cloud with 
a mass well-approximated by the Bonnor-Ebert mass, given 
as (e.g. Palla 2002) 



-1/2 



(23) 



where nfl na i is the final density and T nna i is the final temper- 
ature of the post-shock gas prior to fragmentation. Isobaric 
density evolution gives 

T ps no 



tflnai 



Tflnai 



(24) 



where n is the pre-shock density of the gas and T ps is the 
initial post-shock gas temperature. This initial post-shock 
temperature is roughly 



fcB 



(25) 



where fee is the Boltzmann constant, u sn is the shock veloc- 
ity, and run is the mass of the hydrogen atom. The CMB 
floor defines the final temperature as follows: 



Tflnai = TCMB = 2.7 K(l + z) . 



(26) 



Substituting equations (24), (25), and (26) into equation 
(23) gives 



4M 



fi±fW n ° r 1/2 ( u ^ r 1 (97 ) 

3 V 21 J Vl0 2 cm- : V boOkms-V ' V ' 



Next, as an estimate of a typical preshock density of baryonic 
material in a DM halo, assume the density at the point of 
virialization of the halo, given in Clarke & Bromm (2003) 

as 



no ^ 0.3 cm 



(28) 



Define the characteristic Pop II. 5 stellar mass as 
M char ~ qMbe, with the parameter a — 0.3 describing 
the efficiency with which the mass in the cloud is incor- 
porated into stars. Thus, Mbe gives a robust upper limit 
on the masses of Pop II. 5 stars. This value for a is close to 
that inferred for the formation of stars in the present-day 
Universe, although it is in general a function of the specific 
accretion physics involved in star formation and is not yet 
known with any certainty for the case of primordial star 
formation (e.g. McKee & Tan 2002; Bromm & Loeb 2004). 
Then, using equations (27) and (28), we find 



M cha 



20Mr; 



l + z 
21 



1/2 



(-. 



Wsh 



\200kms- 



(29) 



This formula for the characteristic Pop II. 5 stellar mass is 
the equivalent to equation (9) in Mackey et al. (2003) , which 
was calculated neglecting the effects of HD cooling on the 
shocked primordial gas. Because HD cooling allows shocked 
primordial gas to efficiently cool to the temperature of the 
CMB, the characteristic mass of Pop II. 5 stars is redshift 
dependent, and strong shocks in the primordial gas will thus 
yield higher mass stars, on average, the earlier in the history 
of structure formation they occur. As well, higher velocity 
radiative shocks will tend to yield lower mass Pop II. 5 stars. 
Fig. 9 shows the characteristic Pop II. 5 stellar mass as a 
function of redshift for various SN shock velocities. 
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Figure 9. The characteristic mass of Pop II. 5 stars, formed 
from primordial gas shock-compressed by the first supernovae. As 
shown in equation (29), this characteristic mass varies inversely 
with the velocity of the supernova shock u sn , while increasing 
with redshift. That the masses of these stars explicitly depend on 
redshift is due to the fact that HD cools the shocked primordial 
gas to the temperature of the CMB. For each of the shock veloc- 
ities shown here, the density of the preshock gas was taken to be 
the virial value given in equ. (21). 



4.2 Characteristic stellar mass for DM halo 
collapse 

To calculate the characteristic mass of stars formed in DM 
halo assembly, we assume that the shocked gas, with a post- 
shock temperature T ps roughly equal to the virial tempera- 
ture of the halo, cools isobarically, following the treatment 
of Clarke & Bromm (2003). The pressure of the gas as it 
cools, defined as that necessary to balance the weight of the 
overlying baryons, is given by these authors as 



P vir ~10 7 fc B M 8 2 / 3 (i±i) 4 K 



(30) 



where M x = M/10 x Mq. In taking this value for the pres- 
sure, we have assumed that the mass of the halo is higher 
than the critical mass above which centrifugal effects be- 
come important in slowing the collapse of the gas in the 
halo, given by Clarke & Bromm (2003) as 



M cri t ^ 10 s M ( 



3/2 



An 



(31) 



We here assume that the spin parameter has the fiducial 
value of Ao.i = 1. Again, we have that T nna i = Tcmb = 2.7 K 
(l+z). Thus, the density of the gas when the temperature 
has dropped to T nn ai is given by 



nfinal — Pvir/(fc B TcMB) 



(32) 



Using these values for the final temperature and density in 
equation (23), we obtain the characteristic stellar mass for 
the case of stars formed in collapsing DM haloes as 

1/3, 



M cnar ~ lOMg /o M 



(33) 



where we have again used M cn ar — c*M B e, with a — 0.3. 
The characteristic mass of stars formed in the assembly of 
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DM haloes with masses i; 10 8 M© depends only on the halo 
mass. Thus, stars formed in this process are expected to have 
characteristic masses of the order of 10 Mq, comparable to 
the SN-shocked case. 



5 SUMMARY AND CONCLUSIONS 

We have investigated the importance of HD as a coolant of 
the primordial gas in four distinct cases: the evolution of 
supernova-shocked gas, with shock velocities <; 100 km s _1 ; 
the evolution of gas shocked during DM halo collapse, with 
shock velocities 100 km s _1 ; freely collapsing unshocked, 
un-ionized gas clouds within minihaloes; and freely collaps- 
ing clouds of ionized gas within relic H II regions. In the 
cases where the primordial gas is ionized, HD is an impor- 
tant coolant. Thus, behind strong shocks and in relic H II 
regions, the primordial gas is able to cool by HD line emis- 
sion to temperatures Ss 100 K. In the case of strong shocks, 
HD cooling can lower the temperature of the primordial gas 
to the temperature of the CMB in a fraction of the age of the 
Universe at redshifts z J; 10. In much of the literature, it is 
assumed that metals are required to lower the temperature 
of interstellar gas to the CMB temperature (e.g. Bromm et 
al. 2001; Schneider et al. 2002; Clarke & Bromm 2003). We 
emphasize that through HD cooling gas of purely primor- 
dial composition, without any metals, can cool to the CMB 
temperature, the lowest temperature attainable by radiative 
cooling. All that is required is for a large enough abundance 
of HD to be produced, which we have shown can be done 
under strongly-shocked conditions in which the primordial 
gas is ionized and evolves isobarically. 

We have calculated the characteristic mass of primordial 
stars formed from gas that has cooled to the CMB temper- 
ature floor. We estimate that Pop II. 5 stars, formed from 
gas shocked either through SN explosions or during the as- 
sembly of the first dwarf galaxies, have characteristic masses 
that are redshift dependent and of the order of 10 M©. 

The implications of this ability of the primordial gas to 
effectively cool to low temperatures are numerous. HD cool- 
ing offers a mechanism by which low mass metal-poor stars 
can be formed very early in the history of the Universe, at 
redshifts z <; 10. Thus, HD cooling may play a role in the 
formation of the most metal-poor stars that have been ob- 
served. In particular, it has been suggested that the metal 
abundance patterns of two recently discovered extremely 
metal poor halo stars, HE0107-5240 and HE1327-2326, can 
be explained by these stars having formed from gas enriched 
by the elements released in the type II supernova explo- 
sion of a 20-25 M star (Christlieb et al. 2002; Frebel et 
al. 2005). While the first generation of stars which formed 
in minihaloes of mass 10 6 M© would likely have been far 
more massive than this (<; 100 M©), the characteristic mass 
that we calculate for Pop II. 5 stars falls just in this mass 
range. Therefore, it may have been the supernova explosions 
of Pop II. 5 stars that dispersed metals into the primordial 
gas that later formed the most metal poor stars that arc 
observed today. 

As we have shown, the evolution of the primordial gas 
following the DM halo collapse that formed the first dwarf 
galaxies may have allowed for the formation of primordial 
stars with masses of the order of 10 M©. The formation of 
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Figure 10. The characteristic mass of stars as a function of red- 
shift. Pop III stars, formed from unshocked, un-ionized primordial 
gas are characterized by masses of the order of 100 M©. Pop II 
stars, formed by gas which is enriched with metals, are formed 
at lower redshifts and have characteristic masses of the order of 
1 Mq. Pop II. 5 stars, formed from shocked primordial gas, have 
characteristic masses reflecting the fact that they form from gas 
that has cooled to the temperature of the CMB. Thus, the charac- 
teristic mass of Pop II. 5 stars, given by equation (29), is a function 
of redshift and is of the order of 10 Mq. 



the oldest globular clusters, which formed within ~ 10 yr 
after the big bang, must have required the presence of metal 
coolants to account for the observed abundance of low-mass 
stars in them (see Ashman & Zepf 1992; Kang et al. 1990; 
Bromm & Clarke 2002). 

Overall, the fact that primordial gas can cool to the 
CMB temperature by HD line emission suggests that pri- 
mordial stars with masses at least an order of magnitude 
lower than those characteristic of the first stars could have 
formed. If this is indeed the case, then we obtain a red- 
shift dependence for the characteristic stellar mass of three 
distinct modes of star formation, shown schematically in 
Fig. 10. Pop III stars, formed in the first minihaloes out 
of un-ionized gas that could effectively cool only to ~ 200 
K through H2 line emission, would have been the first stars, 
with masses of the order of ~ 100 Mq. Then, as these first 
stars ended their short lives by exploding as supernovae, 
the strong shocks thus generated would have initiated the 
formation of large fractions of HD, allowing the primordial 
gas to cool to the CMB temperature floor, and so making 
possible the formation of Pop II. 5 stars with masses of the 
order of 10 M© . As well, shocks associated with the forma- 
tion of primordial dwarf galaxies, with masses <; 10 s Mq, 
could have triggered HD cooling, as well, yielding stars with 
masses of the order of 10 M©. Finally, with the dispersion 
of sufficient amounts of metals into the interstellar medium, 
Pop II star formation would have set in, with gas clouds ef- 
fectively cooling to the CMB temperature floor at redshifts 
2^3 and to ~ 10 K up to the present day, the latter tem- 
perature set by dust heating (see Clarke & Bromm 2003; 
Mackey et al. 2003). 

In the present calculations we have neglected radiative 
transfer and the effects of photodissociation on chemical 
abundances. Because the HD photodissociation time is much 
longer than the HD formation time in a LW radiation field 
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produced solely by the shocked gas in the cases we study, we 
expect that it is safe to neglect the effects of photodissociat- 
ing radiation when the only UV flux is from the shocked gas. 
However, with a strong UV background, such as might be 
generated by first generation massive stars in a primordial 
dwarf galaxy, photodissociation of both HD and H2 will be- 
come much more important in determining the thermal and 
chemical evolution of the primordial gas. The effects of HD 
and H2 self-shielding, as well as of HD shielding by H2, will 
have to be carefully included in any calculation of the cool- 
ing of the gas. Thus, more detailed future work will have to 
account for the effects of a possible strong UV background 
and of HD self-shielding, in order to better discern the pro- 
cess of the cooling and fragmentation of shock-compressed 
primordial gas. In particular, if the primordial gas cools in 
an environment with a strong enough UV background, then 
the photodestruction of HD may be so effective as to keep 
Xhd, the abundance of HD, below AnD.crit, and so to pre- 
vent the primordial gas from cooling to the the temperature 
of the CMB within a Hubble time. In this case, metals would 
be required for gas to cool to the CMB temperature floor. 

In addition to a detailed treatment of the photodissoci- 
ation and self-shielding of HD, consideration of the various 
possible heating mechanisms in the high-z Universe, such 
as cosmic rays from the first supernovae, must be included 
in a more thorough treatment of the cooling of the primor- 
dial gas to temperatures approaching the CMB limit (e.g. 
Rollinde, Vangioni & Olive 2005). Our idealized treatment 
serves to illustrate the importance of HD for primordial star 
formation in varying environments, specifically, for the effi- 
cient cooling of shocked primordial gas to the temperature 
of the CMB. To more fully explore the importance of the HD 
coolant, we will carry out three dimensional numerical sim- 
ulations of the fragmentation process in shocked primordial 
gas. 



ACKNOWLEDGMENTS 

We would like to thank Avi Loeb, Paul Shapiro, Simon 
Glover, and Naoki Yoshida for helpful comments and dis- 
cussion. This work was supported in part by NASA Swift 
grant NNG05GH54G 

REFERENCES 

Abel T., Bryan G.L., Norman M.L., 2002, Sci, 295, 93 

Alvarez M., Bromm V., Shapiro P., 2005, ApJ, submitted (astro- 

ph/ 0507684) 
Ashman K.M., Zepf S.E., 1992, ApJ, 384, 50 
Barkana R., Loeb A., 2001, Phys. Rep., 349, 125 
Barsuhn J., 1977, A&AS, 28, 453 

Bond J.R., Arnctt W.D., Carr B.J., 1984, ApJ, 280, 825 
Bromm V., Clarke C.J., 2002, ApJ, 566, LI 
Bromm V., Coppi P.S., Larson R.B., 1999, ApJ, 527, L5 
Bromm V., Coppi P.S., Larson R.B., 2002, ApJ, 564, 23 
Bromm V., Ferrara A., Coppi P.S., Larson R.B., 2001, MNRAS, 
328, 969 

Bromm V., Larson R.B., 2004, ARA&A, 42, 79 

Bromm V., Loeb A., 2004, New Astron., 9, 353 

Bromm V., Yoshida N., Hernquist L., 2003, ApJ, 596, L135 

Christlieb N. et al., 2002, Nat, 419, 904 

Ciardi B., Ferrara A., 2005, Space Sci. Rev., 116, 625 



Clarke C, Bromm V., 2003, MNRAS, 343, 1224 

Draine B. T., 1978, ApJS, 36, 595 

Draine B. T., Bertoldi F., 1996, ApJ, 468, 269 

Flower D.R., Le Bourlot J., Pineau des Forets C, Roueff E., 2000, 

MNRAS, 314, 753 
Frebel A., et al., 2005, Nature, 434, 871 
Galli D., Palla F., 1998, A&A, 335, 403 
Galli D., Palla F., 2002, Planet. Space Sci., 50, 1197 
Glover S. C. O., 2005, Space Sci. Rev., 117, 445 
Haiman Z., Rees M., Loeb A., 1997, ApJ, 476, 458 
Kang H., Shapiro P.R., 1992, ApJ, 386, 432 

Kang H., Shapiro P.R., Fall S.M., Rees M.J., 1990, ApJ, 363, 488 
Kitayama T., Yoshida N., 2005, ApJ, 630, 675 
Kitayama T., Yoshida N., Susa H., Umemura M., 2004, ApJ, 613, 
631 

Larson R.B., 1998, MNRAS, 301, 569 

Le Petit F., Roueff E., Le Bourlot J., 2002, A&A, 390, 369 
Lipovka A., Nunez-Lopez R., Avila-Reese V., 2005, MNRAS, 361, 
850 

Machida M., Tomisaka K., Nakamura F., Fujimoto M., 2005, ApJ, 
622, 39 

Mackey J., Bromm V., Hernquist L., 2003, ApJ, 586, 1 
Mac Low M.-M., Shull J.M., 1986, ApJ, 302, 585 
Madau P., Ferrara A., Rees M.J., 2001, ApJ, 555, 92 
Mathis J. S., 1957, ApJ, 125, 318 
McKee C.F., Tan J.C., 2002, Nat, 416, 59 
Mori M., Ferrara A., Madau P., 2002, ApJ, 571, 40 
Murray S., Lin D., 1990, ApJ, 363, 50 
Nakamura F., Umemura M., 2001, ApJ, 548, 19 
Nakamura F., Umemura M., 2002, ApJ, 569, 549 
Oh S.P., Haiman Z., 2002, ApJ, 569, 558 
Omukai K., Yoshii Y., 2003, ApJ, 599, 746 
Pottasch S. R., 1961, ApJ, 135, 385 

O'Shea B., Abel T., Whalen D., Norman M., 2005, ApJ, 628, L5 
Palla F., 2002, in Maeder A., Meynet G., eds., Physics of Star 

Formation in Galaxies. Springer, Berlin, p. 16 
Rollinde E., Vangioni E., Olive K., 2005, ApJ, 627, 666 
Salvaterra R., Ferrara A., Schneider R., 2004, New Astron., 10, 

113 

Savin D., 2002, ApJ, 566, 599 

Schneider R., Ferrara A., Natarajan P., Omukai K., 2002, ApJ, 
571, 30 

Shapiro P. R., Kang H., 1987, ApJ, 318, 32 
Shull J.M., Silk J., 1979, ApJ, 234, 427 
Storey P.J., Hummer D.G., 1995, MNRAS, 272, 41 
Suchkov A. A., Shchckinov Y. A., Edclman M.A., 1983, Astro- 
physics, 18, 360 
Uehara H., Inutsuka S., 2000, ApJ, 531, L91 
Whalen D., Abel T., Norman M., 2004, ApJ, 610, 14 
Yamada M., Nishi R., 1998, ApJ, 505, 148 

Yoshida N., Abel T., Hernquist L., Sugiyama N., 2003, ApJ, 592, 
645 



© 2005 RAS, MNRAS 000, 000-000 



